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Abstract. Exchange interaction has been studied for electrons in coupled quantum 
dots (QD's) by a configuration interaction method using confinement potentials with 
different profiles. The confinement potential has been parametrized by a two-centre 
power-exponential function, which allows us to investigate various types of QD's 
described by either soft or hard potentials of different range. For the soft (Gaussian) 
confinement potential the exchange energy decreases with increasing interdot distance 
due to the decreasing interdot tunnelling. For the hard (rectangular-like) confinement 
potential we have found a non-monotonic behaviour of the exchange interaction as 
a function of distance between the confinement potential centres. In this case, the 
exchange interaction energy exhibits a pronounced maximum for the confinement 
potential profile which corresponds to the nanostructure composed of the small inner 
QD with a deep potential well embedded in the large outer QD with a shallow potential 
well. This effect results from the strong localization of electrons in the inner QD, which 
leads to the large singlet-triplet splitting. Implications of this finding for quantum logic 
operations have been discussed. 



PACS numbers: 73.21.La 



Exchange interaction in coupled quantum dots 



2 



1. Introduction 

An exchange interaction between electrons localized in coupled quantum dots (QD's) 
is a very promising tool for manipulating qubits in semiconductor nanodevices [JJ. 
This interaction can change the spin of the electron, which allows us to perform the 
quantum logic operations with spin qubits [H El El HI [5]. Recently, the exchange- 
interaction induced spin swap operations in coupled QD's have been simulated by 
a direct solution of a time-dependent Schodinger equation [6]. The quantum logic 
operations can be performed in the laterally coupled QD's [3 [9] and quantum wire 
QD systems pm EH]. Conditions of a realization of the logic operations with qubits 
in QD's are determined by the properties of electron states in QD nanostructures. 
The electron states of the laterally coupled QD's have been studied theoretically in 
Refs. [121 H31 HH [13 HSl UHl HSl Eni EH [22]. In electrostatically gated QD's [7], 
the properties of the electron quantum states can be tuned by changing the external 
voltages applied to the gates [23] . 

In order to perform the high-fidelity quantum logic operations with spin qubits the 
exchange interaction should be possibly strong. The exchange interaction energy, defined 
as the difference between the lowest triplet and singlet energy levels, depends on the 
localization of electrons in the QD's. In general, the stronger the electron localization 
the stronger the exchange coupling. The electron localization is determined by the 
profile of the potential confining the electrons within the QD's. Therefore, the exchange 
interaction depends on the shape and range of the confinement potential. Usually, the 
confinement potential in coupled QD's is modelled by the two-centre parabolic [HI [21] 
or Gaussian [161 EH] potential. 

In the present paper, we propose the two-centre power-exponential (PE) potential 
|24j, which allows us to study a broad class of confinement potentials with different 
shapes. The one-centre PE potential [21] is very well suited for a description of the 
electrostatic QD's [231 [25]. It has a form 

V(r) = -V exp[-tr/R) p ], (1) 

where Vq is the depth of the potential well, R is the range of the potential, which 
determines the QD size, and parameter p > 2 describes "a softness" of the potential, 
i.e., a smoothness of the QD boundaries. If p ~ 2, the potential is "soft", and if p > 4, 
the potential is "hard", i.e., it possesses the walls with large steepness. Parameter p can 
be used to describe a different steepness of the QD boundaries. Therefore, PE potential 
(CE]) can be applied to a modelling of the electrostatically gated QD's [7J [23] and self- 
assembled QD's with compositional modulation [26]. The influence of the smoothed 
interfaces on the electronic and optical properties of GaAs/Al x Gai_ x As QD's have been 
studied in a recent paper [27] . 

The exchange interaction can be controlled by internal nanostructure parameters, 
e.g., size and geometry of the coupled QD's, and external electric and magnetic fields 
[TBI El [20]. It has been shown [16] that the asymmetry of QD's leads to the increase of 
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this interaction. Recently, the size effects in the exchange coupling have been studied 
for quantum wire QD systems 

In the present paper, we investigate the influence of the shape and range of the 
confinement potential on the exchange interaction in laterally coupled QD's using the 
two-centre PE potential. The paper is organized as follows: a theoretical model is 
presented in Section II, the calculation methods and results are given in Section III, 
Section IV includes the discussion and Section V - the conclusions and summary. 



2. Theory 

We study the system of two electrons confined in laterally coupled QD's with identical 
confinement potentials. The geometry of the nanostructure, which consists of two 
separated QD's, is illustrated in Figure 1. We describe each QD by the two-dimensional 
(2D) rotationally symmetric potential well centred at positions ±a, where a = (d/2,0) 
and d is the distance between the centres of the potential wells. 




Figure 1. Schematic of the coupled QD system. Circles A and B represent the QD's, 
d is the distance between the centres of the confinement potentials, a is the position 
vector of the centre of the right QD. 



In the effective mass approximation the Hamiltonian of the system reads 

H = h 1 + h 2 + , (2) 

4irE E s r 12 

where hj (j = 1,2) is the one-electron Hamiltonian, Eq is the electric permittivity 
of the vacuum, e s is the static relative electric permittivity, r\ 2 = \r\ — r 2 \ is the 
electron-electron distance, and = (xj,yj) denote the electron position vectors. The 
one-electron Hamiltonian has the form 
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where m e is the electron effective band mass and V(vj) is the confinement potential. 
We assume that the effective mass and the static electric permittivity do not change 
across the QD boundaries. This assumption is well justified for the electrostatic QD's 
based on GaAs [23]. For the coupled QD's the confinement potential is taken on as a 
sum of PE potentials [Eq. (CQ)] centred at ±a. In the explicit form, 

V(r) = -V { ex V [-(r A /RY} + exp[-(r B / R) p }} , (4) 

where r& = |r + a| and t~b = |r — a|. 

Formula (jl]) defines a broad class of two-centre confinement potentials with different 
shapes. For fixed potential well depth V and range R parameter p characterizes the 
softness (hardness) of the confinement potential. The smaller (larger) p the potential 
is more soft (hard). For p = 2 the confinement potential has the Gaussian shape. This 
potential is soft, i.e., the potential walls at the QD boundaries, have fairly small steepness 
and are partly penetrable for the electron. For p > 4 the confinement potential becomes 
hard, i.e., the potential walls at the QD boundaries are steep. For p > 10 we deal with 
the rectangular-type, very hard, confinement potential. Figure 2 shows the different 
profiles of the two-centre confinement potential [Eq. (jl])]. In Figure 2 and throughout the 
present paper we are using the donor Rydberg R D = m e e A / \?>2ii 2 E 2 ] E 2 s h 2 ) as the unit of 
energy and the donor Bohr radius an = ^ireoh 2 / (m e e 2 ) as the unit od length. For GaAs, 
Rd ~ 6 meV and ~ 10 nm. We note the essential difference in potential profiles 
between the soft (p = 2) and hard (p = 10) confinement potentials for intermediate 
distances d, i.e., for d ~ 3 ar, in Figure 2. In the case of soft confinement potential, 
the resulting two-centre potential changes smoothly with increasing d from the single 
to the double potential well. In this case, the increasing intercentre distance d leads 
first to the increase of the potential well size and next to the formation of the potential 
barrier for d > 3 an- If the confinement potential is hard, then - for intermediate 
intercentre distances d - we obtain a narrow deep potential well surrounded by fairly 
wide potential steps, on which the potential is flat [cf. Figures 2(b,c)]. This shape of 
the confinement potential corresponds to the core-shell QD-nanostructure, which was 
realized in CdTe/CdS self-assembled QD's [28]. In this case, we deal with the compound 
QD nanostructure, which consists of the small inner QD embedded in the larger outer 
QD. 

3. Results 

We are mainly interested in the influence of the shape of the confinement potential, in 
particular its softness, on the electronic properties of coupled QD's. Therefore, we have 
performed the calculations for fixed depth V = 25 Rp and range R = 2 a D (except 
otherwise specified). These values of the parameters correspond to laterally coupled 
GaAs QD's [9]. 
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Figure 2. Confinement potential profile for different distances d between the potential 
well centres as a function of coordinate x. Solid curves correspond to the hard 
confinement potential (j> = 10), dashed curves correspond to the soft confinement 
potential (p = 2). The plots for Vq = 25B,£> and R = 2ajj. 



3.1. One- electron problem 

The one-electron eigenvalue problem for the one-centre confinement potential has been 
solved for the Gaussian potential in Ref. [25] and for the PE potential in Ref. [24|. 
In the present paper, we consider the one-electron bound states for the two-centre 
confinement potential (j3J). We have solved the eigenvalue problem for Hamiltonian 
([3]) by the imaginary time step method [30] applying the finite-difference approximation 
of Hamiltonian ([3]) on the 2D grid with 101 x 101 mesh. In this case, the accuracy of 
the method [30] is high enough to treat the obtained numerical solutions as exact. 

Figure 3 shows the results for the six lowest-energy levels, which correspond to 
the one-electron states used to a construction of two-electron configurations. Figures 
3(a) and 3(b) display the results for the soft (p = 2) and hard (p = 10) confinement 
potential, respectively. For d = we deal with the single QD with the potential well of 
double depth (2Vb), while for large d the QD's are separated by the potential barrier. 
For intermediate d the superposition of hard-wall potentials corresponds to the inner- 
outer QD nanostructure [cf. Figures 2(b) and 2(c)]. These properties of the confinement 
potentials affect the one-electron states (Figure 3). The one-electron energy levels are 
monotonically increasing functions of interdot distance d. For d = the excited-state 
energy levels exhibit a degeneracy, which is removed for d > 0. We note that for 
intermediate values of d, i.e., for d ~ 3 ap, the energy levels quickly increase with 
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Figure 3. Lowest-energy levels Eo, ■ ■ . , E5 of the single electron confined in coupled 
QD's as functions of distance d between the confinement potential centres for (a) soft 
confinement potential with p = 2 and (b) hard confinement potential with p = 10. 

The dependence of the ground-state energy Eq on the softness parameter p is 
depicted in Figure 4. If the confinement potential is more hard, i.e., for large p, the 
ground-state energy takes on the lower values, which means that the electron is more 
strongly bound due to the larger quantum "capacity" of the QD. 

3.2. Two-electron problem 

Two electrons confined in the coupled QD's form molecular-like states, called the 
artificial molecules p2J [13]. We solve the two-electron eigenvalue problem by the 
configuration interaction (CI) method using the one-electron numerical solutions 
obtained in the previous subsection. Augmenting the calculated spatial wave functions 
by the eigenf unctions of the z component of the electron spin we obtain the one- 




Figure 4. Ground-state energy Eq of the single electron confined in the coupled QD's 
as a function of interdot distance d for different softness p of the confinement potential. 



electron spin-orbitals ip ua (r), where v and a are the orbital and spin quantum numbers, 
respectively. Next, we construct Slater determinants 

X„Oi,r 2 ) = ^[-0 nCT1 (ri)Vv 2CT2 (r 2 )] , (5) 

where A is the antisymmetrization operator and n labels different two-electron 
configurations with well-defined total spin. According to the CI method the two-electron 
wave function is a linear combination of Slater determinants ([5]) 

#(r x , r 2 ) = c„XnOi, r 2 ) . (6) 

n 

In the present calculations, we have used 15 Slater determinants, which were constructed 
from the six lowest-energy one-electron states. We have checked that including 20 Slater 
determinants in expansion ([6]) only slightly improves the results, but the computation 
time considerably increases. Two-electron Hamiltonian (T5]) has been diagonalized in 
basis ([6]). All the matrix elements, including the electron-electron interaction energy, 
have been calculated numerically on the 2D grid defined for the one-electron problem. 

The results for the lowest-energy levels of the singlet (E$) and triplet (E T ) states 
are displayed in Figure 5, which also shows the exchange interaction energy Ex energy 
defined as 

Ex = E T — E s . (7) 

The energy of the singlet and triplet states increases monotonically with increasing 
distance d between the confinement potential centres. Therefore, the binding energy, 
defined as E B = —E T; g, of both the states is a decreasing function of interdot separation, 
i.e., the two-electron states are more weakly bound for more separated QD's. Figure 
5 shows the different behaviour of the the exchange energy for the soft and hard 
confinement potential. If the confinement potential is soft [cf. Figure 5(a) for p = 2], 
the exchange energy decreases with increasing intercentre distance and for large d is 
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Figure 5. Lowest-energy levels of the singlet (Es) and triplet (Et) states (left scale) 
and exchange energy Ex (right scale) as functions of distance d between the centres 
of the confinement potential for (a) soft (p = 2) and (b) hard (p — 10) confinement 
potential. 

negligibly small. This behaviour results from the decreasing interdot tunnelling with 
increasing distance between the QD centres. If the confinement potential is hard [cf. 
Figure 5(b) for p — 10], the exchange energy increases with intercentre distance (for 
small d), exhibits a pronounced maximum for d ~ 3.3aD (for p = 10), and next rapidly 
falls down to zero for large d. We shall discuss the origin of this behaviour in the next 
subsection. 

Figure 6 shows the dependence of the exchange energy on intercentre distance d 
and softness p of the confinement potential. We see that the maximum of the exchange 
energy is more pronounced if parameter p is large, i.e., the confinement potential is 
hard. Moreover, the maximal value of the exchange energy is larger if the confinement 
potential is more hard. 
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Figure 6. Exchange energy Ex as a function of distance d between the confinement 
potential centres for different values of softness parameter p. 



3.3. Electron density distribution 

In order to explain the different behaviour of the exchange energy for the soft and hard 
confinement potentials, we have calculated the one-electron probability density 



and displayed it in Figures 7, 8, and 9. 

The properties of the electron density are determined by the shape of the total 
confinement potential, which is different in the x and y directions. The ^-dependence of 
the one-electron probability density (Figure 7) is very similar for the singlet and triplet 
states. In the x direction, the electrons are localized in the single potential well for 
small d (including d — 0). For large interdot separation (cf. Figure 7 for d — 4.2od) 
the electrons are localized with equal probabilities in the two quantum wells separated 
by the energy barrier. At the intermediate distances, i.e., for d ~ 3ao, the confinement 
potential profile corresponds to the core-shell QD nanostructure with the inner and 
outer QD's (cf. Figure 7 for d = 2.8a£> and 3. 2a d). In this case, the electrons are 
mainly localized within the inner QD with the deep potential well. Figure 7 shows that 
- for all interdot distances - the electron localization in the singlet state is stronger then 
in the triplet state. In particular, we observe that - for d — 3.2ar> - the triplet electron 
density is fairly large in the outer QD region. 

Figure 8 displays the shape of the confinement potential and the one-electron 
probability density as functions of y for x — 0. For d < 3ap we observe the essential 
qualitative difference in the y-dependent electron distribution between the singlet and 
triplet states. In the singlet state, the electron distribution in the y direction is similar to 
that in the x direction. However, in the triplet state, the electron density distribution is 
more spread out and exhibits two pronounced maxima separated by a minimum at y — 0. 
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Figure 7. One-electron probability density p\ (solid curves) and confinement potential 
V (dashed curves) for p = 4 as functions of x for y = for the singlet (left panel) and 
triplet (right panel) states. 



For sufficiently large interdot distances (d > 3a D ), the electron density distributions for 
the singlet and triplet states are again very similar. 

The comparison of Figures 7 and 8 for d < 3ao allows us to extract the information 
about the distribution of the electrons in the triplet state. The electrons with the same 
spin are strongly localized in the inner QD in the x direction (Figure 7), but are more 
weakly localized in the y direction (Figure 8). The weaker localization of the electrons in 
the y direction results from the larger effective range of the potential well in this direction 
(cf. dashed curves in Figures 7 and 8). Figures 7 and 8 show the reasons of the existence 
of the maximum exchange interaction for the hard confinement potential. We see that 
the increase of the exchange coupling is caused by the strong electron localization in 
the inner QD and the considerable difference of the electron density distributions for 
the singlet and triplet states. These effects lead to the increase of the triplet-singlet 
energy difference, which in turn results in the formation of the pronounced maximum 
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Figure 8. One-electron probability density p\ (solid curves) and confinement potential 
V (dashed curves) for p = 4 as functions of y for x = for the singlet (left panel) and 
triplet (right panel) states. 



of exchange energy at intermediate distances between the confinement potential centres 
(cf. Figure 5). 

These effects are additionally illustrated in Figure 9, which displays the contour 
plots of the one electron density distribution and confinement potential on the x — y 
plane. The arrows schematically show the sites with the largest electron density for 
the singlet and triplet states. The plots for d = and d = 2.8a£> show that - in the 
triplet state - each electron is localized at a different site. The plots for d = 2.8ao 
correspond to the maximum of the exchange energy for p = 4. The triplet electron 
density distribution is anisotropic on the x — y plane. For small d the electrons with 
the same spin are aligned in the y direction. This configuration rapidly changes for 
d ~ 3a D and the triplet density distribution becomes more spread out in the x direction 
(cf. Figure 9 for d = 3.2od). Figures 7, 8, and 9 show that - for large interdot distances 
- the electron density distribution is the same in the singlet and triplet which leads to 
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Figure 9. Contours of the one-electron density distribution for the singlet (left panel) 
and triplet (middle panel) states, and the corresponding confinement potential profiles 
(right panel) for different distances d between the potential well centres in the case of 
the hard confinement potential with p = 4 and R = 2ajj. Full (empty) arrows in the 
right panel show the sites with the maximum electron density for the singlet (triplet) 
states. 

the singlet-triplet degeneracy, i.e., the exchange energy tends to zero for large d. 

Figure 10 displays the dependence of maximum exchange energy E^ ax (cf. Figure 
6) on parameter p and range R of the confinement potential. The maximum exchange 
energy increases with increasing parameter p, i.e., increasing hardness of the confinement 
potential, and is the largest for the rectangular-like potential well. However, the 
maximum triplet-singlet splitting quickly decreases with increasing range R, i.e., 
increasing QD size. The independence of E^ ax allows us to determine the size effect in 
the exchange interaction [IT] . This dependence can be parametrized as follows: 

Ex ax {R) =a/R + b, 

with a = 9.975 and b = —0.073 (in donor units). 

For fixed parameters p and V$ the maximum exchange energy depends not only on 
R, but also on distance d between the centres of the confinement potential. Therefore, 
the intercentre distance d = d max , which corresponds to the maximum of the exchange 
energy, changes along the curves E^ ax (R) and E^ ax (p) in Figure 10. The dependence of 
d max on parameters R and p is plotted in Figure 11. It appears that intercentre distance 
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Figure 10. Maximum exchange energy E^ ax as a function of range R of the 
confinement potential (solid curve for p = 4) and softness parameter p (dashed curve 
for R — 2an) calculated for Vq — 25 Rd- 

d max is a linear function of confinement potential range R. Figures 10 and 11 allow us 
to determine the parameters of the QD nanostructure, for which the exchange energy is 
maximal. As we have pointed out these parameters correspond to the inner-outer QD 
nanostructure (cf. Figures 2 and 6). For the coupled QD's separated by the potential 
barrier the exchange energy is considerably smaller. 

4. Discussion 

The parametrization of the confinement potential given in Eq. (jlj) is sufficiently flexible 
to model various types of QD's, among them the electrostatically gated QD's [7] and 
self-assembled QD's with the compositional modulation [26]. Using the confinement 
potential (JH) we can describe the effects of smoothness of the QD boundaries [27J. 
For the intermediate separations d between the potential- well centres (d ~ 3ci£> for 
R ~ 2a d) the resulting potential profile corresponds to the core-shell QD nanostructure 
with the attractive core potential well. This profile of the confinement potential is 
characteristic for the compound QD nanostructure with the inner and outer QD's, 
which has been recently realized by the chemical growth of self- assembled CdTe/CdS 
QD's [28]. The inner-outer QD nanostructure can also be realized in the laterally 
coupled electrostatically gated QD's [7J. The corresponding profile of the confinement 
potential can be obtained by applying the suitably tuned external voltages to the two 
pairs of different gates. When varying the external voltages, we can tune the shape of 
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Figure 11. Distance d max between the confinement potential centres, for which the 
exchange energy is maximal, as a function of confinement potential range R (solid 
curve for p = 4) and softness parameter p (dashed curve for R = 2a_o). 

the confinement potential, which in turn leads to a desirable change of the exchange 
interaction. 

We have shown that the exchange energy is maximal for the compound QD's, which 
consists of the inner QD with the deep potential well and outer QD with the shallow 
potential well. The maximum of the exchange energy is caused by the strong electron 
localization in the inner QD and the large difference of the electron localization between 
the singlet and triplet states (cf. Figures 7, 8, and 9). This result gives us a possibility of 
designing the nanostructures, in which the exchange interaction is sufficiently large to be 
used in quantum logic gates [6]. The strong exchange interaction is very important for 
performing high-fidelity quantum logic operations with spin qubits in QD's [TJ El [3j [6] . 
The recent study [6] of the exchange-interaction induced spin swap operation shows 
that the swapping the electron spins is very effective, i.e., the electron spins are fully 
interchanged in the possibly short time, if the confinement potential changes from the 
double potential wells separated by the barrier to the deep potential well located inside 
the shallow potential well. The later potential profile corresponds to the inner-outer 
QD nanostructure [cf. Figure 2(b,c)]. 

The exchange interaction in the coupled QD's can also be tuned by applying the 
external magnetic [H3, [T7J EH] and electric [8] fields. The increasing magnetic field causes 
the decrease of the exchange interaction energy since it lowers the energy of the triplet 
state and leads to the singlet-triplet degeneracy for high magnetic fields [HJ [17] . The 
rapid changes of the external electric field enabled Petta et al. [8] to perform a coherent 
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manipulation of spin qubits in laterally coupled QD's. Szafran et al. [15] showed that 
the asymmetry of the QD's considerably increases the exchange energy. The size effects 
in the exchange coupling were studied by Zhang et al. The increasing size of the 

QD nanostructure leads to the decrease of the exchange interaction due to the decreasing 
interdot tunnelling [llj. Using the present results we can also determine the size effect 
in the exchange coupling. In particular, Figure 10 shows that the exchange energy scales 
as ~ 1/R with the increasing size R of the QD's. 

5. Conclusions and Summary 

In the present paper, we have studied the lowest-energy singlet and triplet two-electron 
states in laterally coupled QD's and determined the exchange interaction between the 
electrons. The application of the two-centre PE function parametrization [Eq. (j4j)] 
enabled us to investigate a large class of realistic confinement potentials with different 
shapes. We focus on the dependence of the exchange energy on the distance between 
the confinement potential centres and also on the shape and range of this potential. The 
dependence of the exchange energy on intercentre distance d is qualitatively different for 
the soft (p = 2) and hard (p > 4) confinement potentials. For p = 2 the exchange energy 
is a monotonically decreasing function of d, while for p > 4 the exchange energy increases 
with d for small d, reaches the maximum for intermediate d, and next decreases to zero 
for large interdot separations. This knowledge allows us to predict the nanostructure 
parameters, in particular their size and geometry, which maximize the exchange energy. 

In summary, we have found that the exchange energy is maximal for the confinement 
potential, which corresponds to the compound QD nanostructure, consisting of the inner 
QD with deep potential well embedded in the outer QD with shallow potential well. 
The corresponding core-shell confinement potential can be obtained in the form of the 
inner-outer QD nanostructure realized in self-assembled QD's and electrostatically gated 
QD's. We have also investigated the tuning of the exchange interaction by changing the 
parameters of the coupled QD nanostructure and pointed out the importance of the 
present study for the quantum logic operations with electron spins. 
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